Thermal Operator Representation of Finite Temperature Graphs II 
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Using the mixed space representation, we extend our earlier analysis to the case of Dirac and 
gauge fields and show that in the absence of a chemical potential, the finite temperature Feynman 
diagrams can be related to the corresponding zero temperature graphs through a thermal operator. 
At non-zero chemical potential we show explicitly in the case of the fermion self-energy that such 
a factorization is violated because of the presence of a singular contact term. Such a temperature 
dependent term which arises only at finite density and has a quadratic mass singularity cannot be 
related, through a regular thermal operator, to the fermion self-energy at zero temperature which is 
infrared finite. Furthermore, we show that the thermal radiative corrections at finite density have 
a screening effect for the chemical potential leading to a finite renormalization of the potential. 

PACS numbers: ll.10.Wx 



I. INTRODUCTION 

In an earlier paper (referred to as I), we gave a 
simple derivation of an interesting relation be- 
tween finite temperature Feynman graphs and the cor- 
responding zero temperature graphs within the context 
of a scalar field theory in real as well as in imaginary 
time formalisms. We showed that the derivation is par- 
ticularly simple if one uses a mixed space representation 
of the graphs in the (t,p ) space [H ul HI an d the proof 
of the correspondence is particularly simple in the closed 
time path formalism 0, 0, [ljj. Explicitly, for any TV- 
point graph (at any loop) in a scalar field theory at finite 
temperature, the relation can be written as 
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with E{ = Vk 2 + m 2 , rii — n(Ei) denoting the Bose- 
Einstein distribution function associated with the inter- 
nal propagators and Si — S(Ei) represents a reflection 
operator that changes Ei — > —E\. In / characterizes 
the number of internal propagators, V the total number 
of vertices in the graph (with the usual relation for the 
number of loops L = I — V + 1) and 5^(k,p) enforces 
the conservation of momentum at the vertex v. We de- 
note the internal and the external three momenta of a 
graph generically by k,p respectively. Furthermore, 7^ 
represents the integrand of the finite temperature graph 
(after the internal time coordinates have been integrated 



in the mixed space) so that it has the dependence 

7 l P(TX,t a ), (3) 
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where t a ,a = 1,2, . . . ,N denote the external time coor- 
dinates of the graph while 7^ is the integrand of the 
same graph at zero temperature (with the internal time 
coordinates integrated). The operator relating the 
integrands of the two graphs was termed the thermal op- 
erator and the most important property of this operator 
is that it is independent of time coordinates and carries 
the entire temperature dependence of the (finite temper- 
ature) graph. This interesting result is calculationally 
quite useful and allows us to study directly many ques- 
tions of interest at finite temperature such as Ward iden- 
tities and analyticity [ill Il2t Il3| . In I, we had shown 
that this simple relation arises as a consequence of the 
factorization of the finite temperature propagator in the 
scalar field theory into a basic thermal operator acting 
on the zero temperature propagator and had studied var- 
ious properties associated with this thermal operator. In 
particular, we had shown that the basic thermal operator 
for the propagator corresponds to a projection operator 
that projects onto the space of periodic functions. (We 
recall that while the finite temperature propagator for the 
scalar field satisfies periodic conditions following from the 
Kubo-Martin-Schwinger condition |14j . the zero temper- 
ature propagator does not.) For a complex scalar field 
with a chemical potential, on the other hand, we showed 
that the basic thermal operator is much more complex in- 
volving time derivative terms. In this case, we could not 
give a general proof of a thermal operator representation 
such as although we showed, for specific complicated 
graphs, that a nontrivial factorization nonetheless arises. 

In this paper, we extend our analysis in I to theories 
involving fermions as well as gauge theories. The analy- 
sis for gauge theories is particularly of interest since the 
interaction terms (non-Abelian three point interaction as 
well as the interaction of the ghost fields) involve deriva- 
tive terms. We find in all cases that if there is no chemical 
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potential present, a thermal operator representation for 
finite temperature graphs naturally follows. On the other 
hand, for a fermion with a chemical potential, as in the 
case of the complex scalar field discussed in I, the basic 
thermal operator is complicated involving time deriva- 
tives and we find that a thermal operator representation 
for graphs fails. This failure is traced to the fact that in 
such theories, the self-energy develops a quadratic mass 
singularity because of radiative corrections at finite den- 
sity. The paper is organized as follows. In section II, wc 
discuss fermion theories at finite temperature (without a 
chemical potential) and show that the basic factorization 
of the thermal propagator arises much as in the scalar 
field theory. The proof of the thermal operator repre- 
sentation for an interacting theory involving scalar and 
fermion fields is direct in the closed time path formalism 
which we discuss. In section III, this analysis is extended 
to gauge theories where we show that the basic factoriza- 
tion of the thermal propagator leads to a thermal opera- 
tor representation for any graph at finite temperature in 
spite of interaction terms involving derivatives. The ther- 
mal operator representation is explicitly worked out for 
the contribution of the ghost loop to the self-energy of the 
gauge field. In section IV, we study an interacting the- 
ory of gauge fields and fermions with a chemical potential 
(for example, QED at finite density) and show that in this 
case the basic factorization of the thermal propagator for 
the fermion involves a dependence on time derivatives. In 
this case, the basic thermal operator can also be written 
equivalently as one without a time derivative but with a 
matrix structure. We work out the fermion self-energy 
in this theory explicitly and show that a thermal opera- 



tor representation fails. We trace this failure in section 
V to the fact that the quantum corrections in this the- 
ory lead to a quadratic mass singularity at finite density. 
By analyzing the pole of the fermion propagator in this 
theory, we show that the chemical potential has a finite 
renormalization due to radiative corrections and we dis- 
cuss some interes ting aspects of this phenomenon (see, 
for example, ref. |l5| for a discussion from the point of 
view of the renormalization group evolution). We con- 
clude with a brief summary in section VI. In appendix 
A, we study the 0+1 dimensional Chern-Simons QED to 
bring out some interesting features of the thermal opera- 
tor representation in lower dimensions while appendix B 
describes briefly the derivation of some of the formulae 
used in the text. 



II. FERMIONS (WITHOUT CHEMICAL 
POTENTIAL) 

In this section, we will study an interacting theory of 
scalar and fermion fields at finite temperature (without 
chemical potential). As we have shown in I, the proof 
of the thermal operator representation for any graph is 
direct in the closed time path formalism jRl^llOjl . There- 
fore, for simplicity, we will discuss the theory in this 
formalism although everything we say also holds in the 
imaginary time formalism. Indeed, in reference 0] the 
validity of the thermal operator representation has been 
shown to hold when there is no chemical potential. Let us 
consider the theory described by the Lagrangian density 
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The factorization of the scalar propagator has already 
been discussed in I and we simply recapitulate here the 
essential results. In the closed time path formalism, the 
propagator has a 2 x 2 structure which can be written as 



A^(t,E) ^O^A^^E), 
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where the 2x2 matrix structure of the propagator is 
labelled at any temperature as 



\A {T l(t,E) A iT l(t,E) J 
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and the basic thermal operator is the scalar operator 

O b T) (E) = l + n B (E)(l-S(E)), (7) 



where 



n B (E) = 



CT — 1 



, E = E(M) = y / p 2 + M 2 . (8) 



The components of the propagator at zero temperature 
have the following explicit forms in the mixed space 
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2E 



where the operator L(e) takes the limit e — > 0. 

For fermions, on the other hand, we know that the 
components of the 2x2 matrix propagator (in the closed 
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time path formalism) at finite temperature have the mo- 
mentum space representation 



(10) 



with (L(e) is the operator taking the limit e — > intro- 
duced earlier) 
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S?l(p) = 2n{j> + m) (e(-p Q ) - n F (\p \)) S(p 2 - m 2 ), 
S^(p) = 2-K{j> + m) (9(p ) - n F (\ PQ \)) S(p 2 - m 2 ), 
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Here n.F ( |po | ) represents the Fermi-Dirac distribution 
function 
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and the temperature dependent terms reflect the anti- 
periodicity condition satisfied by the fermion propagator. 

The components of the fermion propagator in (flip can 
be Fourier transformed in the energy variable to give 



S^(t,p) = 
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where E = E(m) = \/p 2 + 



and 



OP(E) = l-n F (E)(l-S(E)). 
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The components of the zero temperature propagator have 
the explicit forms 

— \e{t)A{E) e -^ E - %e)t 
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where 
A(E) = j°E-j-p+m, 



B{E) = -j E-j-p+m. (16) 



It is worth remarking here that, as in the case of the 
scalar propagator, it is easy to verify that the basic ther- 
mal operator in l|14|) is a projection operator, namely, 



(17) 



and in the present case projects onto functions satisfying 
anti-periodicity properties. 



Thus, we see that in spite of a matrix structure (from 
the Dirac gamma matrices) of the fermion propagator, 
the thermal propagator factorizes in terms of a basic 
thermal operator (|14fl which is a scalar quantity much 
like in the case of the scalar field theory. Furthermore, 
it is independent of the time coordinates and as a re- 
sult, the thermal operator representation for any graph 
can be obtained as follows. First, let us suppose that 
we have a graph with only external vertices (that is, a 
one loop graph). A typical N- vertices graph involving 
fermion propagators will have the general form shown in 
Fig. (The external vertices can be of "±" type, but we 
choose all of them to be of "+" for illustrative purposes 
only. The same derivation will go through for vertices 
of any type since the basic thermal operators in (JJJ and 
(|14f) are scalar quantities and have the same form for 
any component of the propagator.) In this case, at finite 
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FIG. 1: A typical one-loop diagram involving fermions (solid 
lines) and scalar fields (dashed lines). For simplicity, the ver- 
tices are all assumed to be of "+" type. 



temperature, the value of the graph can be written as 
(we set g — 1 = A and ignore all the multiplicative fac- 
tors coming from the vertices for simplicity. The external 
momenta are all assumed to be flowing into a vertex, and 
the internal momentum fej flows from vertex i to vertex 
z + l and we identify ijv+i = ii, kisr+i = fci,pjv+i = Pi-) 
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where we have identified the thermal operator for the 
graph as 
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The basic thermal operators (E) , F ; (E) are de- 
fined in Eqs. J7J and l|14|l respectively and we have to 
remember that in 1)21(1 



Ei = 



k 2 + m 2 for i = 1,2, . . . , J F 



(22) 



'k 2 + M 2 for £ = I F + 1,...,JV 
As a result, we can write 

fi ( 2 -) 3 

showing that in this case, the finite temperature graph 
can be given a thermal operator representation. Further- 
more, since in the closed time path formalism the range of 
time integration at finite temperature continues to be the 
same as at zero temperature and since the basic thermal 
operators J7J and l|14f> are independent of time coordi- 
nates (so that they can be taken outside the integral), 
such a factorization of any graph with internal time co- 
ordinate (that need to be integrated over) continues to 
hold and in general, for any iV-point graph (at any loop) 
at finite temperature, we have the thermal operator rep- 
resentation 

r - T) -/ IlTSlW^MO^, (24) 

1=1 ^ ' V — l 

where the thermal operator follows from Q21J to be 

0^=f[0^( El ) f[ 0g\ Ei ), (25) 

i=l i=J F + l 



with Jp j I representing respectively the number of inter- 
nal fermion propagators and the total number of internal 
propagators. 



III. GAUGE THEORIES 

We have seen thus far that the thermal operator rep- 
resentation for any Feynman graph at finite tempera- 
ture holds for theories involving scalar and fermion fields 
(without a chemical potential). However, physically 
gauge theories are more interesting and in this section 
we will discuss a non-Abelian gauge theory at finite tem- 
perature. Let us consider a Yang-Mills theory (where 
the gauge fields belong to SU(n)) in the Feynman gauge 
described by the Lagrangian density 
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£ = -^ F %E^ a - - (d ■ A a y + dTDpf, (26) 

where a = 1,2,.. . ,n 2 — 1 and (we set the coupling to 
unity for simplicity) 

D»c a = d^c a + f abc Alc c , 
F; v = dpAZ - d„A; + f abc AlA c v . (27) 

In this case, in the closed time path formalism, the gauge 
and the ghost propagators at finite temperature have the 
momentum space representation 

KTiP) = *,f3 = ±, (28) 

(T) 

where A a p(p) represent the components of a massless 
scalar propagator at finite temperature and have the ex- 
plicit forms 



a2(p) 

A {T l(p) 
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-L(e). 
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+ 2irn B (\p Q \)6(p 2 ) 



where L(e) is the limiting operator introduced earlier. 

By taking the Fourier transform of 1|29|) with respect 
to po, we can obtain the components of the gauge and 
the ghost propagators in the mixed space. We already 
know from © and that at finite temperature the com- 
ponents of the scalar propagator factorize in the mixed 
space representation. It follows, therefore, that the com- 
ponents of the gauge and the ghost propagators also fac- 
torize in the mixed space representation as 
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where the same basic thermal operator Og (E) leading 
to a factorization of the gauge and the ghost propagators 
coincides with that for a scalar propagator defined in 
(with E = \p | for a massless field). Furthermore, all the 
components of the propagator factorize in the same man- 
ner and the basic thermal operator is independent of the 
time coordinate. It is worth remarking at this point that 
we have chosen to work in the Feynman gauge for simplic- 
ity. In any other covariant gauge fixing, only the Lorentz 
structure of the gauge propagator generalizes, but the ba- 
sic factorization of the thermal propagator continues to 
hold. Furthermore, we consider the case of a vanishing 
chemical potential here for simplicity. 

Since the gauge and the ghost propagators factorize in 
the same way as in a scalar field theory, the thermal op- 
erator representation of any graph at finite temperature 
would seem obvious. However, unlike in a scalar field 
theory, the interactions in a non-Abelian gauge theory 
involve time derivatives in the mixed space (for exam- 
ple, the three gluon vertex or the ghost interaction ver- 
tex) and, in principle, may complicate the general proof 
of the thermal operator representation of an arbitrary 
graph. On the other hand, we note that the basic ther- 
mal operator (JTJ in the factorization of the propagators 
(|3U|I is independent of time coordinates. As a result, it 
follows trivially that 

d t OP{E) = Og\E)d t , (31) 



so that the basic thermal operators in a propagator can 
be trivially commuted past the derivatives in the vertices 
leading to a thermal operator representation of any ar- 
bitrary graph. Let us illustrate this with the example 
of the one loop graph (see Fig. |5J) depicting the ghost 
contribution to the gauge self-energy. In this case, the 
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FIG. 2: The ghost loop contribution to the self-energy of the 
gauge boson. The "dots" represent the action of the time 
derivatives. 



graph can be written as (we suppress all multiplicative 
factors associated with the diagram for simplicity, iden- 
tify &3 = k\ and look only at the contributions coming 
from a time derivative in the vertex which may present 
a challenge to a thermal operator representation) 
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where in the last step we have identified the thermal op- the thermal operator representation holds simply be- 

erator for the graph to be cause the basic thermal operator commutes with the time 

m(T) /n( T ) r \mi T )i derivative operator. With this observation, it is clear that 

— b \ U*/b \ 2 >- for an interacting non-Abelian gauge theory, we can write 

Thus, we see that in spite of the presence of time deriva- the thermal operator representation for any arbitrary N- 

tives in the vertices in the mixed space representation, point graph at finite temperature (involving gauge and 
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ghost vertices of "±" type) at any loop as 

r; ,.. QN( T) = [f[** f[(2n)%V(k, P hP^ 

= I ii(|ii( 2i )^^) o,T, ^"" (o ' ( 34 ) 

i— 1 ^ ^ v — 1 

where we have suppressed the Lorentz indices associated 
with the graph and the thermal operator for the graph 
has the form 

0^=f[0^(E i ). (35) 
i=i 

In a similar manner, it can be shown that in the absence 
of chemical potentials any diagram in an interacting the- 
ory involving gauge fields, fermions and scalar fields at 
finite temperature will have a thermal operator represen- 
tation. The interesting and challenging case, however, 
seems to be in the presence of a chemical potential which 
we discuss in the next section. 



IV. FERMIONS WITH A CHEMICAL 
POTENTIAL 

Let us next consider QED at finite temperature and 
density. In the Feynman gauge, the theory is described 



by the Lagrangian density 



C = -\f^F^ + i> (ilp- m) ^-\{d- A) 2 + m ^tV, 

where /i represents the chemical potential associated with 
the fermion and the covariant derivative is defined to be 



= - ieA^. (37) 



In H36fl we have neglected the free Lagrangian density for 
the ghosts which is not relevant for our discussions. 

As we have argued before, there is no chemical poten- 
tial associated with the photon. As a result, the prop- 
agator for the gauge boson in the Feynman gauge will 
continue to be what we have discussed in the last section 
(without any internal indices). On the other hand, in 
momentum space in the closed time path formalism, the 
components of the fermion propagator in the presence of 
a chemical potential take the forms 



S {T f } 

o( t >a0 
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(38) 



where E = \/p 2 + m 2 . Equation Ij38(l clearly reduces to 
when = 0. The Fourier transform of (|38|l in the 
Po variable leads to the propagator in the mixed space 
which can be seen to have a nontrivial factorization (as 
is the case with the complex scalar field discussed in ref. 
I) 



(39) 



S^\t,p) = e^O^\E)Sff\t,p), 

where, as before, a, /3 = ± and the basic thermal operator 
in l|39(l has the form 



\E) = 1- 



-(l-S(E)) 



+ s m ** t (40) 



where we have defined 



and S^g°\t,p) represent the components of the fermion 
propagator at zero temperature and zero chemical poten- 
tial given in l|15|) . 

We note that the structure of the basic thermal opera- 
tor in 14()p is quite analogous to the case of the complex 
scalar field with a chemical potential discussed in ref. I. 
This is a scalar operator which involves a time derivative 
operator (it does not depend on the time coordinate). 
Unlike in the case of the scalar field, however, in this 
case, we can equivalently define a basic thermal operator 
which is independent of the time derivative, but instead 



n F (E ± (j,), 



(41) 
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is a matrix (in the Dirac space), namely, in this case we 
can also write 



Here A(E), B(E) denote the Euclidean rotation of the 
matrices in <|16ll . namely, 
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where 
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°a/3 



= 1 - (n^A(E) - n+B(E)) - S(E)), (43) 



(t,i>), (42) A{E) = ij E-f -p+m, B(E) = -ij Q E-^-p+m. (48) 
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where A(E), B{E) are matrices defined in 1)16(1. The two 
forms of the basic thermal operator are related through 
the first order equation satisfied by the fermion propa- 
gator in the mixed space. However, the scalar form of 
the basic thermal operator (in spite of the time deriva- 
tive operator) is easier to use than the matrix one and we 
will carry out our discussions in terms of the factoriza- 
tion H39|) ■ We note that both forms of the basic thermal 
operator in l|4U|l and 143( 1 can be checked to be projection 
operators, namely, 



(7» 



(44) 



and enforce the necessary anti-periodicity in the present 



Let us next analyze the one loop fermion self-energy 
graph at finite temperature and density in QED (see La- 
grangian density 1(36(1 ) at one loop in the imaginary time 
formalism. The two point function in Fig. [3] can be ex- 
plicitly evaluated and leads to (we note that the fermion 
self-energy is simply the two point function with momen- 
tum conserving delta functions factored out and we iden- 
tify &3 = fci in the derivation below) 




P x fe v P, 

FIG. 3: The fermion self-energy at one-loop order. 



As in the case of the complex scalar field with a chem- 
ical potential (in ref. I), here we note that the presence 
of the time derivative term in the basic thermal operator 
(|40|l makes it difficult to give a general proof of a thermal 
operator representation for any arbitrary graph, partic- 
ularly in cases involving internal time coordinates which 
are integrated over. In the scalar theory, we had shown in 
a specific complicated graph that a (nontrivial) thermal 
operator representation results in spite of the presence of 
the time derivative operator and this led to the hope that 
a thermal operator representation may result in general. 
In the fermionic theory, we will show through an explicit 
calculation that the thermal operator representation fails 
in the presence of a chemical potential and we will trace 
this failure (in the next section) to a renormalization of 
the chemical potential because of radiative corrections. 
For the purpose of explicit calculations, we will use the 
imaginary time formalism |5j, llo . Il7| where there is no 
doubling of fields. In this case, the factorization of the 
thermal propagator (|39|l can be written as 

S^( T ,p) = e» T O^\E)S^' \T,p), (45) 

where O 1 ^'^ (E) is the rotation of the operator in H4U|) 
to imaginary time 

0^\E) = 1-!£±^f (1-5(25)) 

_!tZk(l+5(f5))|, (46) 

and the propagator at zero temperature in the absence 
of a chemical potential has the form 
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S (m (r,p) = — [e(r)A(E)e-^ + 9(-r)B(E)e E ^] 
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d3k ~(Tm) 
(2tt) 3 72 ' 
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where we have identified the integrand with j^'^ (to 
avoid confusion with the Dirac gamma matrices) and 
have carried out the integration over k 2 using one of the 
delta functions (and have identified k\ — k). The inte- 
grand has the explicit form 



72 



-e 2 ^S^{T - t'M )i^ t \t - t 1 , Ei) 
: -e 2 \e^ T - T '^0 { p{E 1 )6 { p^ } (E 2 ) 

x ( 7 ^ (0 ' 0) (T-T',fc 2 )^D^(t-t',E 1 )) 



2i(n+(E 2 )-np(E 2 )) l0 



6(t - t') 



E( Ei 



2ie 2 {n+(E 2 )~n ¥ (E 2 )) l0 



El -El 



S(t-t'), (50) 



where the modified fermion operator O^'^ is similar to 
the one in l|46|l except that it involves d T j (E\ + E 2 ) and 



(47) 



E x = \k\, E 2 = \j{k+p) 2 +m 2 , pi=p, k 2 = k+p. 

(51) 

Thus, we see explicitly from (|5U|I that, unlike the case 
of the complex scalar field discussed in ref. p|, here a 
thermal operator representation for the fermion two point 
function at one loop breaks down in the presence of a 
chemical potential. The additional term leading to the 
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breakdown of the thermal operator representation is a 
contact term which vanishes for u — > since 



lim (n+(E 2 )-n F (E 2 )) 



0. 



(52) 



Such an identity has already proved quite useful in the 
solution of the + 1 dimensional Chcrn-Simons QED 8] 
and with the explicit form of the fermion propagator in 
(|45|l it can be checked that this is true. It also follows 
from that 



V. RENORMALIZATION OF CHEMICAL 
POTENTIAL 

Let us note that the chemical potential can be thought 
of as a constant background electrostatic potential. This 
is particularly clear if we note that the free part of the 
fermion Lagrangian density in l|36|) can be written in the 
Euclidean space as 



(53) 



As a result of this structure of the theory, one can de- 
rive certain identities as follows. Adding the sources for 
the fields, we can write the generating functional for the 
theory in (|36[1 (in imaginary time) as (we set h = 1) 

Z[ f M,J fl ,r,,fj]=e- w ^ J ^ 



This, in turn, leads to the identity 



dW 
da 



dr I d 3 x I / 



SW SW 



8n(x) ^° 8n(x) 
S 2 W 



Srj(x)Srj(x) 



(55) 



By taking the second order derivative with respect to 77, fj 
and setting all sources to zero, this leads to the identity 
(when we Fourier transform the spatial coordinates) 

dS^(ri-r 2 ,p) 
dji 

= i\ T dr S^(n - r,p )7 S (3>) (t - r 2 ,p ).(56) 
Jo 

I 



1 gns<T,n) 
nl dfi n 



= {irS^ 7 o---S^hoS^\ (57) 



where there are n insertions of 70 on the right hand side 
and we have used a compact notation suppressing the 
internal time integrations. 

Using the above result, we can see that a correction 
to the chemical potential may arise from successive in- 
sertions of the operator Su 70 in the fermion propagator 
,£( T >W. The resulting corrected propagator S r ( T >/ i +** t ) is 
then obtained by summing the geometric series 



Y,Zo(i) n S^ (^70) ---S^ (Su l0 ) = 

^ Wafts™ (TM+Su) 

^ W nl du n 

n=0 ^ 



(58) 



Thus, the effect of such insertions is to shift the chemical 
potential to an effective value given by p + Sp,. In order 
to evaluate this shift (finite renormalization) we have to 
perform a more systematic analysis of the fermion self- 
energy. To this end, let us calculate the complete self- 
energy in momentum space [BJ (namely, Fourier trans- 
form l|49(l in the time variable and factor out overall 
energy-momentum conserving delta functions) 



dr e lpoT ^ T '^(T,p). 



(59) 



A direct evaluation yields 



d 3 k 



(2tt) 3 4£i£ 2 



-B 



A 



(1 



(l + n B 0Ei)-ME 2 )) (MEi) + n F {E 2 )) 



E 1 +E 2 -(ip + 
n B {E 1 )-n+{E 2 )) 



p) Ei- 
(n B (Si) 



E 2 + (ip + p) 



E l +E 2 + (ip + p) E l -E 2 - (ip + p) 



(60) 



where we have identified 

A = 2(i"/ E 2 - 7 • k 2 - 2m) 

B = 2(-i~f E 2 -j-k 2 -2m). (61) 



It is clear from the explicit form of Q60[) that the self- 
energy (in fact, any thermodynamic function) is a func- 
tion of (ipo + p, p) where the extra p dependence comes 
from the explicit dependence of the fermion distribution 
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functions on the chemical potential. The analytic contin- 
uation of the self-energy can be carried out appropriately 
[iH through 

iPo=*Po+M->-R> = (Po + AO(l+*e), ( 62 ) 
say for the time ordered amplitude. We recognize from 



(|60|l that the real part of the coefficient of the 7 term in 
(|6U[) coincides exactly with the coefficient of the contact 
term in the mixed space in (|5t)|l when po + /j = 0. 

In fact, let us note that after analytic continuation, we 
can write the coefficient of the 7 term in (|60(l as 



1 n it,a, . e f d 3 k 1 rr 

-Tr 7 °£( I » (p) = / - — — <^ 



2 J (2tt) 3 E x 
+n B (Pi) 



1 



E 1 +E 2 -P E 1 +E 2 + Pq 
1 



Ei + E 2 - P Ei+E 2 + P E1-E2- P E 1 -E 2 + P Q 
n+{E 2 ) + n F (E 2 ) ( 1 1 1 



1 



2 \E 1 +E 2 -P E 1 +E 2 + P Q E 1 -E 2 -P E 1 - E 2 + P Q 

n+{E 2 ) - n F (E 2 ) f 1 1 1 1 

E l +E 2 - P E l +E 2 + P E x -E 2 - P E x - E 2 + P 



(63) 



Let us denote the terms involving the braces as Si and 
the terms in the last parenthesis as E 2 (after integration) 
so that we can write 

l Tr7 o E (7» (p) = EiCPo,^) + £ 2 (P 0) p»- (64) 

From the explicit form of the terms in (|63|) . we note 
that under Pq — » — Po,/i — > — /U, both £i,£ 2 change 
sign. However, since £1 is manifestly symmetric under 
/.t — > it is anti-symmetric under Po — > — Po and, con- 
sequently, vanishes when P = 0. On the other hand, £ 2 
is anti-symmetric under fi — > —fi and, consequently, is 
symmetric under Pq — > — Po and does not vanish when 
Po = and, in fact, yields the contact term in (|50H . Using 
the results in appendix B, we can evaluate £ 2 explicitly 
in the high T limit and when p = 0, it has the form (for 
T»m) 



£ 2 (Po = 0,p = 0) 



eV(7r 2 r 2 + M 

67T 2 TO 2 



(65) 



This has a quadratic mass singularity that arises only 
when fi =/= (and, therefore, at finite charge density). 
The self-energy at zero temperature and chemical poten- 
tial on the other hand, is infrared finite (logarithmic 
infrared divergences only appear if one expands around 
the singular point p 2 — to 2 ). Consequently, the singu- 
lar term l|65|) cannot be related to the zero temperature 
fermion self-energy through a regular thermal operator. 
We believe that the presence of this strongly divergent 
infrared behavior is responsible for the failure of the ther- 
mal operator representation in the case of a nonvanishing 
chemical potential. 

To understand the renormalization of the chemical po- 
tential, we have to analyze the poles of the fermion prop- 
agator. We note that with the self-energy corrections, the 



complete two point function at one loop (in Minkowski 
space) can be written as 

liS^^y 1 ^) = i> - m + M7° - E(p, T, /x), (66) 

where £ can only be calculated in some limit such as the 
high T limit. The analysis of the poles of the propagator 
in this limit, even in the absence of a chemical potential, 
is highly nontrivial. It is known in the absence of a chemi- 
cal potential ^3] that in the leading order at high temper- 
ature (T 3> to), the fermion propagator has an absolute 
pole at p = 0,po — ±m/ where to/ = represents the 
thermal mass of the fermion. For nonzero p, the fermion 
propagator has only partial poles corresponding to two 
quasi-particle modes. Here we will follow the same anal- 
ysis restricting ourselves to only the absolute pole in the 
leading order at high temperature. Also, we will disre- 
gard those thermal corrections which yield a finite renor- 
malization of the vacuum fermion mass, because such 
non-leading terms are not relevant for the analysis of the 
renormalization of the chemical potential. In the leading 
order at high temperature (T 3> to), the terms in l|6U|) 
which are even under fi — > — \i lead to 



"»/7°. Po 



In- 



P 



2p P -p 

P V 2p P 



P 



(67) 



where we have defined p — \p\ and p denotes the unit 
vector along p. In the presence of a chemical potential 
we have 



2 e T 2 , l L 
m f = T ( T + ^2 



(68) 
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which is well known (see, for example, [17]'), The non- 
analytic behavior of E e ven at Pq = 0,p = is obvious 
from Hfi7[) . We also note that all the 7 • p terms in Ij67(l 
vanish in the limit p = (basically because in this case, 
there is no direction available to contract the gamma 
matrix) so that we can write 



(f=0) ~ -/ 7 ° = 7 °Si(Po,P = 0). 



P 



(69) 



The absolute pole in the fermion propagator continues to 
be at p = and at this point, the terms proportional to 
70 in l|60(l which are odd under (j, — > —fx yield 



odd 



e 2 (P 2 +m 2 )7° 
2^ 

= 7°£ 2 (P ,P = 0). 



dfcfc 2 



(n?(fc)-n+(*0) 



4fc 2 P 2 



(P 2 - m 2 ) 2 
(70) 



Using the results in appendix B, this can be evaluated 
in the high temperature limit and shows that at |Po| = 
\M\ ~ eT m it is well behaved and has the value 



odd 



(P = M,p = 0) = 7 Q £ 2 (P = M,p = 0) 



e 2 /^7° 
8tt 2 ■ 



(71) 



With these results, the analysis of the pole when p = 
becomes quite straightforward. We note from l|66|) that 
when p — 0, the propagator will have a pole provided 



7°P) 



>(P = 0) 



,)dd 



(p = 0) = 0. (72) 



We note from (|68|l and l|69|) that at very high temperature 
mj > m so that the fermion mass may be neglected in 
the above equation. All the other terms are proportional 
to 7 . If we expand £ dd around Pq = M, the equation 
(|72|l takes the form 



Pn 



m 



P> 



e 2 /x 
8^ 



(P -M)£ odd (P = M)- 



0. 



(73) 



Here we have used (|71|l and £ odd denotes the derivative 
of S odd with respect to Pq. The root of this equation 
and, therefore, the location of the pole is given by 



e 2 u 

P = M = T7 ^ ± m f 



1 



e 2 /i 2 



16tt 2 3 \ 327r 2 (vT 2 T 2 
which can be equivalently written as 



Po 



9 

e \x 

167T 2 



±m/ 1 



327r 2 (vr 2 T 2 + /i 2 ) 



(74) 



(75) 



To the order that we are working, this can be simplified 
to give the location of the pole at 



or. 



Po + AM 1 



Pa + A'R = ± 



2\/2 



±m f 



(76) 



T 2 + a. 



Here we have identified 

MR = M [ 1 - 



16tt 2 



(77) 



which can be interpreted as the renormalized chemical 
potential due to the radiative corrections of the theory. 
Since it is associated with a physical pole of the propaga- 
tor, we expect this result to be gauge independent which 
we have explicitly checked. Such a finite renormalization 
of the chemical potential has the effect of screening the 
chemical potential because of thermal interactions at a 
finite density. This is consistent with our earlier obser- 
vation that the chemical potential can be thought of as 
a constant Abelian electrostatic potential and Abelian 
gauge fields lead to a screening effect. 



VI. SUMMARY 

In this paper, we have extended our analysis of the 
thermal operator representation for Feynman graphs at 
finite temperature to theories involving fermions as well 
gauge fields. We have shown that as long as there is 
no chemical potential, the thermal operator representa- 
tion holds. We have also discussed in an appendix how a 
thermal operator representation naturally arises in + 1 
dimensional Chern-Simons QED. However, in QED at 
finite temperature and density (nonzero chemical poten- 
tial), we have shown that such a factorization is violated 
because of the appearance of singular contact terms. This 
is explicitly worked out in the case of the fermion self- 
energy at one loop. The reason for this failure of the 
thermal operator representation is traced to the presence 
of a quadratically divergent thermal infrared singularity 
in the self-energy for a non-zero chemical potential (fi- 
nite density). In this case, we find that the chemical 
potential undergoes a finite renormalization due to ra- 
diative corrections. The renormalized chemical potential 
is determined from an analysis of the pole of the fermion 
propagator at high temperature and shows that the ra- 
diative corrections lead to a screening of the chemical 
potential. This is argued to be consistent with the ob- 
servation that a chemical potential can be thought of as 
a constant electrostatic potential and screening is a phe- 
nomenon associated with Abelian gauge fields. 

In conclusion, we would like to point out that the lack 
of a complete factorization in the presence of a chemical 
potential may be related to our choice of generalizing the 
basic thermal operator in terms of the simple reflection 
operator S(E). Finding an alternate basic thermal op- 
erator possibly dependent on other nontrivial operators 
and determining its consequences on factorization is an 
interesting issue which is presently under study. 
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APPENDIX A: 0+1 DIMENSIONAL 
CHERN-SIMONS QED 

As we have shown in ref. I as well as in this paper, a 
thermal operator representation for any finite tempera- 
ture graph holds in any theory in the absence of chemical 
potential. The reflection operator S(E) in the thermal 
operator simply changes E — > — E and in this way incor- 
porates the negative energy contributions into the graph. 
Furthermore, the thermal operator is independent of the 
time coordinate which plays a significant role in the gen- 
eral proof of the thermal operator representation. All of 
this is true in higher dimensional field theories. How- 
ever, in a 0+1 dimensional field theory (quantum me- 
chanics), the situation is different because the energy is 
positive and the question arises as to whether a thermal 
operator representation holds for such a theory as well. 
Furthermore, as is well known, in Chern-Simons QED in 
0+1 dimensions amplitudes beyond the one point func- 
tion vanish at zero temperature while all the higher point 
functions are nonzero at finite temperature yj- There- 
fore, it is interesting to analyze how the nonzero finite 
temperature amplitudes in such a theory can arise from 
a thermal operator acting on trivial amplitudes. 

Let us recall that the Lagrangian for the 0+1 dimen- 
sional Chern-Simons QED in the Euclidean space is given 
by 

C = <4>{d T -iA + m)i/j -inA, (Al) 

where n represents the Chern-Simons coefficient and we 
have set the coupling for the gauge field to unity for sim- 
plicity. In this case, the fcrmion propagator (in the imag- 
inary time formalism) has the form 

S (T ' m) (T) = e-' mT (0(t) -n F (m)) 

= e- mT O { p' m) (T)S^ \T), (A2) 

where we have identified the basic thermal operator of 
the theory with 

Of m) (r) = l-n F H(l + S(r)). (A3) 

There are two things to note here. First, in the 0+1 di- 
mensional case, the mass term corresponds to a chemical 
potential and second, the basic thermal operator con- 
tains a reflection operator that reflects the time coordi- 
nate and, therefore, is manifestly time dependent. This is 
quite different from the higher dimensional cases we have 
studied where the basic thermal operator is independent 
of the time coordinate. In the 0+1 dimensional theory, 
on the other hand, we do not have higher loop diagrams 
(photon is non dynamical) and, consequently, the time 
dependence of the basic thermal operator does not pose 




FIG. 4: A typical iV-point graph in + 1 dimensional CS 
QED. 

a problem in deriving a thermal operator representation 
for any graph. 

For any graph of the form in Fig. 0] we can immedi- 
ately write down the thermal operator representation as 
(using the identification rjy+i = t% and the fact that the 
exponential term around a closed loop vanishes) 

r-\N N 

r£' m) = -§rIP (T ' m) (^ + i) 

i=i 

N 
i=l 

= G^r%°\ (A4) 
where we have identified 

N 

0^ = l[0^ m \n-n +1 ). (A5) 
i=i 

Indeed we see that formally there is a thermal opera- 
tor representation for any graph in the 0+1 dimensional 
Chern-Simons QED. 

In practice, on the other hand, we know that in this 
theory 

r^ 0) = 0, for N > 2. (A6) 

The way (|A4|I works in practice is as follows. Let us 
identify 

= Ti — Ti + 1, (A7) 

and define the amplitude in (IA4fl in the limiting manner 

r^ m) = Km o£' m) r% 0) . (A8) 

The limit tjv+i — > T\ is assumed to be taken only at the 
end of the calculation (after the action of the thermal 
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operator). With this, let us show explicitly how the cor- 
rect finite temperature two point function arises from a 
trivial zero temperature amplitude. 

We note that the two point amplitude at zero temper- 
ature is given by 

p (o,o) = _^ 5 (o,o) (tli9)5 (o,o) (t2i3) = l 9 (r 13 )0fo t3 ). 

(A9) 

This, of course, vanishes if we identify T3 = t% (or 
T 2,3 = t~2.i). However, we are not supposed to take the 
limit until the thermal operator has acted on the zero 
temperature amplitude. Letting the thermal operator 
act on (IA9I). we obtain 



i(r,m) 



f' m) (ri2)Of' m, (T 2 3)kri2)fi(T 2 3) 



lim O 



lim - (9{t 12 ) - n F (m)) (6>(r 23 ) - n F (m)) 



= -~n F (m) (1 -n F (m)). 



(A10) 



This is indeed the correct finite temperature result (for a 
single fermion flavor) [|| and this shows how the thermal 
operator correctly reproduces the nonzero finite temper- 
ature amplitudes from trivial zero temperature ones if 
the operation is carried out in a limiting manner. 



APPENDIX B: DERIVATION OF THE HIGH T 
LIMIT 

In this appendix, we evaluate some integrals which are 
used in the text. Let us first consider 



/>oo 

^oddO ti > m > r ) = / dk (n F (E - n) - n F (E 
Jo 



■M)) 



(Bl) 

which is essential for obtaining the high temperature 

limit in i|71|) . Here E = y/k 2 + m 2 and in the high tem- 
perature limit T 3> m, it is possible to obtain a closed 
form expression for for an arbitrary /1 in the follow- 
ing way. Let us expand the integrand in a power series in 
fi and integrate term by term. Every term in this series 
is well behaved and leads to 



T 



0) = 2/z^ 



fi 2i d 2e n F (t) 



1=0 



{2t+l)\ dt 2i 



-O 



(B2) 

Since the distribution function for the fermion can be 
expanded as [lflj 



n F (t) 



1 1 ^ , s t n 

n=0 



(B3) 



where E(x) represents the Euler polynomials, the expres- 
sion (|B2() can be evaluated in terms of the Euler func- 
tions. A further simplification results from the fact that 



£?o(0) = l, £W(0) =0for^>0. 



(B4) 



Consequently, in this limit, i|B2(l has the form 

n&(*£-Q) = * + o(£). (B8) 

This result is true for any value of \x. 

However, when ^ is also small, one can determine the 

next order correction to S7^ d as follows. We recall an 
alternative expansion of the fermion distribution function 
as 



1 1 



e 2 + l 2 ^ (2n + 1) 2 tt 2 + z 2 ' 



(B6) 



If we substitute this expansion into (|B1J| and regularize 
the inte gral by multiplying p~ e with e — > taken at the 
end |20l l2l|. the integrand can be expanded in a power 
series in S, In this case, each term in the series can 
be integrated and the series can be summed to give Rie- 
mann 's zeta function C(2n+1, 5) 19]. A straight forward 
calculation leads to 



OO 



M ^(-l)«C(2n + l,-)iJ„_ 1 (r 2 )(— ) , 



n=l 



(B7) 



where r 2 = ^ and H n {r 2 ) are polynomials of order n in 
r 2 . For example, for the first few we have 



H (r 2 ) = l, ff 1 (r 2 ) = ~(3 + 4r 2 ),. 



(B8) 



In the limit, ^ — > 0, we recover l|B5|l 

Using similar techniques, we can furthermore show 
that 



dkk 2 (n F (E - fi) - n F (E + /i)) 



|(. 2 T 2 + M 2 )+0(^), (B9) 



0^ n (/i,m,T) = / dkk (n F (E - fi) + n F (E + fi)) 
Jo 

1 /tt 2 T 2 ,\ m x 
= 2(— +^J+°(t)- ( B1 °) 

One may extend the set of formulae (|B5|) . 1B9|I and 
(|B10j) with the help of the basic integral: 



i(p, y) 



^r^ydx = -T(p + l)IA v+1 (-e u ), 

(Bll) 



which can be obtained by expanding the integrand in 
powers of e u and then integrating term by term. Li n (z) 
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is the polylogarithm function, which is the analytic con- 
tinuation to the whole complex z plane of the series (valid 
for n > 1 and \z\ < 1) 



Li n (z) = 



Z K 



k=l 



Using a generalization of the method employed by Haber 
and Weldon in appendix A of reference |2l|, one can 
find a power series expansion of the function Li p+ i(— e"), 
leading to 



where £(z) is the Riemann zeta function, and where the 
singular numerator corresponding to n — p must be in- 
terpreted as its limiting value, In 2. The formula is valid 
for all real values of p > — 1. Notice that for integer p the 
series in (|B12j) . starting from the power u p+2 , contains 
powers of the same parity only, because the zeta function 
vanishes at all negative even integers. 



/(p, u ) = r(p+i)]T 



(i-a^Kfr + i-n)^ (B12) 
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